YITP-00-45 



Gravitational waves in cosmological models of 
Hofava-Witten theory 



Osamu Seto and Hideo Kodama 
Yukawa Institute for Theoretical Physics, Kyoto University 
Kyoto 606-8502, Japan 
(February 1, 2008) 



Abstract 



We study the behavior of gravitational waves and their backreaction on the 
background in cosmological solutions of the five-dimensional Hofava-Witten 
theory. As a dynamical background, we consider two cosmological solutions 
with spatially flat expanding FRW branes, called (j)- and (j)-solutions, in 
which the orbifold size increases and decreases in time, respectively. For these 
background solutions, the wave equation for the tensor perturbation can be 
solved by the method of separation of variables, and the mode functions are 
classified by a separation constant which can be regarded as a graviton mass. 
We show that the spatial behavior of the mode functions are the same for both 
background solutions, but the temporal behavior is significantly different. 
We further show that for the (j) solution, the background bulk geometry 
is unstable against the backreaction of the perturbation, while for the {{) 
solution, the backreaction on the bulk geometry can be neglected. We also 
show that, in contrast to the effect to the bulk geometry, the backreaction 
of the perturbation significantly alters the intrinsic geometry of the brane for 
the (|) solution. 
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I. INTRODUCTION 



Over the past few years, a considerable number of studies have been made on the brane- 
world scenario in which our universe is realized as a boundary of a higher dimensional 
spacetime fllHl6|l. 111 particular, inspired by the recent progress in the heterotic M-theory 
0, five- dimensional brane- world models in which 3-branes are embedded in an effective 
five-dimensional spacetime compactified on S' 1 /Z 2 j3],H have attracted much attention. For 
example, cosmological solutions in the five- dimensional Hofava-Witten theory have been 
discovered |5|-|9[ . The five-dimensional models of Randall and Sundrum, which were proposed 
to solve the Hierarchy problem [15] and to demonstrate an alternative to compactification 
fl3|| , also have many similarities to the five-dimensional Hofava-Witten theory. 

In the brane-world scenario, all ordinary matter fields are confined on the brane, while 
graviton can propagate in the fifth-dimension. Hence, in order to test the idea of the 
brane-world, one needs to study the nature of the gravity in this scenario. The behavior of 
linearized gravity in the Randall - Sundrum models has recently been studied in detail. It 
has been shown that massless modes of the metric perturbation are decoupled from massive 
modes and the Einstein gravity is recovered at low energy scales (See, e.g. [|12] -|i5|]). However 
these investigations have been done only in the highly symmetric background models such 
that four-dimensional maximally symmetric branes, i.e., Minkowski brane [0,0] and de 
Sitter brane ||15|| , are embedded in a five- dimensional anti-de Sitter spacetime, which also is 



locally maximally symmetric in five- dimensions. 

On the other hand, many people have discussed the possibility of a homogeneous and 
isotropic Friedmann- Robert son- Walker (FRW) brane-world in the Hofava-Witten theory 
0-0] and in Randall and Sundrum's scenario [jnj. However, there have so far been few 
attempts to investigate perturbations on such a dynamical cosmological brane-world, al- 
though the four-dimensional real universe has inhomogeneous fluctuations as is shown by 
the observations of the CMB anisotropy. 

As pointed out in Ref . |T7| , it is expected that all gravitons on the brane become massive 



in a dynamical brane-world model. Then, the excitation of massive graviton would become 
a crucial defect in the brane-world scenario, or provide a new model of dark matter in the 
brane-world cosmology. It is therefore important to study the evolution of perturbations on a 
dynamical brane model in order to explore the cosmological consequences of the brane-world 
idea. 

Recently, formalisms for cosmological perturbations on brane-world have been developed 
by several authors [I7|-25j]. In particular, for general cosmological brane-world models, the 
evolution equations for the metric and matter perturbations in the bulk and the boundary 
conditions for them at the brane have been established by Kodama et al. [T7] in terms of 
gauge-invariant variables. The perturbations on a brane are inevitably coupled to the per- 
turbations on the bulk. The evolution equations and boundary conditions for cosmological 
perturbations, in particular, the scalar and vector perturbations on the cosmological brane 
are too complicated to solve. On the other hand, as shown in | JL7|| , when the anisotropic 
stress perturbation vanishes, the tensor perturbations decouple from the matter perturba- 
tion and the boundary condition becomes a Neumann type. Hence, as far as the tensor 
perturbations are concerned, the problem is easier to deal with. 

In this paper, as a first step to investigate the perturbations on a dynamical brane-world, 
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we study the behavior of tensor perturbations on two cosmological brane-world models in 
five-dimensional Hofava-Witten theory. In general, the equations for perturbations are no 
longer separable for a dynamical 3-brane which is not maximally symmetric as a hyper- 
surface. Fortunately, for the cosmological solutions in the Hofava-Witten model found in 
IIU, the evolution equation of tensor perturbations becomes separable. However the evo- 
lution equation is still rather complicated to solve exactly. So, in this paper, we consider 
the late time behavior of the tensor perturbation by using WKB approximation to analyze 
the evolution equation. Further we discuss the backreaction problem and the stability of 
the background solutions as well as the brane motion using the second-order perturbation 
theory. We shall show that one of the cosmological solution is unstable if the backreaction 
of the tensor perturbation is taken into account. We shall also show that the backreaction 
from massive modes significantly alters the brane motion, i.e., the evolution of the brane 
universe, for the solution whose bulk geometry is stable against the tensor perturbation. 

The present paper is organized as follows. In the next section we briefly recapitulate 
the five-dimensional Hofava-Witten theory and the cosmological solutions which are used as 
the background for perturbation. In Sec. Ill we give the equations of motion for the tensor 
perturbation and the boundary condition for them, and analyze the behavior of the tensor 
perturbation by using WKB method. Then in Sec. IV we discuss the backreaction of the 
perturbation on the background. Section V is devote to conclusion and discussion. 

II. COSMOLOGICAL SOLUTIONS OF HORAVA-WITTEN THEORY 

Hofava and Witten have shown that the strongly coupled limit of the E$ x Eg heterotic 
string theory has been identified with M-theory compactified on a S l /Z2 orbifold with Eg 
gauge fields on each orbifold fixed plane 0. After a compactification on a Calabi-Yau 
three-fold, the fields of the standard model can be confined to 3-brane ||. Thus one has an 
effective model in which our four- dimensional universe is a 3-brane embedded in an effective 
five-dimensional spacetime compactified on S' 1 /Z 2 . As was shown by Lukas, Ovrut, Stelle 
and Waldram M , the bosonic sector of this effective model contains the 5-dimensional metric 
9mn-i a modulus V describing the variation of the Calabi-Yau volume along the orbifold, and 
the U(l) gauge field Am and two charged scalars (a, £) parametrizing the anti-symmetric 
tensor field. If we assume that the gauge field and the two charged scalar fields vanish, the 
five-dimensional effective action for the bosonic sector of the Hofava-Witten theory is given 
by i: 




where K5 is the five- dimensional gravitational constant, a is a constant, and is the five 
dimensional spacetime bounded by the branes and M.\ . 

The cosmological solutions with flat FRW branes for this effective action have been 
constructed by Lukas et. al. |J and been generalized to include closed and open FRW branes 
by Reall |J. In the present paper, as a simple case, we shall consider only the model in 
which expanding flat FRW branes are embedded in the five-dimensional bulk. 
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Let x 4 = y be a coordinate in the orbifold direction with y G [— irp, np] and Z2 acting on 
S 1 by y — > —y. The orbifold fixed planes are located at y = 0, 7rp. Then, starting from the 
effective five- dimensional action (|I]) with the metric ansatz: 

ds 2 = a 2 e 2A{t) C(y)(-dt 2 + S^dx'd^) + e 2B{t) D(y)dy 2 , (2) 

one obtains the cosmological solutions of five-dimensional Hofava-Witten theory with flat 
FRW branes §: 

C{y)=D(y)W = ^a\y\ + l, (3) 

e 2A = ^-S^ e 2B = t 2 Sj (4) 

where ao are constants and 5 = ±\/3/2. The field V for these solutions is given by 

V = e B ^C{yf. (5) 

Hereafter we shall refer to the upper and lower choices of sign as the (|) and (J,) solutions, 
respectively. 

The (|) solution describes the model that the four-dimensional FRW universe expands 
while the orbifold space shrinks. On the other hand, the (|) solution describes the model 
in which both the four-dimensional FRW universe and the orbifold space expand, and the 
latter expands faster than the former does. 

For a while, we shall focus our attention on the four-dimensional brane at y = 0. In 
terms of the cosmic proper time r defined by 

e Adt = J??" t <s-M». (6) 

o 



t = a 

The Hubble parameter is given by 



H(t) = = l^lr^/ 2 . (7) 
a dr 2a 

Then the wave-number kn whose wavelength corresponds to the horizon scale is given by 

k H = aH= . (8) 

The scale factor of the four- dimensional FRW brane is written as 

[3-5 

a(r) = a ^_ r| . (9) 

Comparing this scale factor with that of the no-extra-dimension cosmology, we find that the 
brane expands as if it were a standard four-dimensional flat FRW universe which contains 
a perfect fluid obeying the equation of state p = wp with 

^W^y < 10 > 

although the brane- world solutions considered here are vacuum solutions. In this picture, the 
FRW brane in the (|) solution looks like an unphysical universe because the dominant energy 
condition violates, w > 1 . On the other hand, in the (!) case, the FRW brane describes a 
physical universe in the sense that the dominant energy condition holds, < w ~ 0.37 < 1. 
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III. BEHAVIOR OF TENSOR PERTURBATIONS 



A. Wave equation for tensor perturbations 

Gravitational perturbations in the bulk are decomposed into components of the three 
types, scalar, vector and tensor, and each component can be expanded by tensor harmonics 
of the same type on the 3-space of constant curvature. Then, the tensor perturbations 
represent gravitational wave modes in the four-dimensional FRW brane. 

For the action ([I]), the anisotropic stress perturbation vanish both in the bulk and in 
the brane. Then, for the background metric Eq.(||]) and (§), the equation of motion for the 
tensor perturbation and the boundary condition are given, respectively, by 

, . • n • 2 2(A(t)-B(t)) 

H T + (2A(t) + B(tj) H T - + k 2 H T = 0, (11) 

H T ' = 0, aty = 0,7rp, (12) 



where Ht is the expansion coefficient of the metric perturbation, Sg^ = 2a^e 2A C{i))HTT. 



in terms of the tensor-type harmonic tensor Ty on the flat 3-space |17]. Here, —k 2 is an 



eigenvalue of the d'Alembertian on the flat 3-space, and dots and primes denote derivatives 
with respect to t and y, respectively. 

Note that the boundary condition flI2p is simply written as the derivative with respect 
to y as in the static brane case. However, the evolution equation ( pTTf ) contains an additional 
friction term BH?, which does not exist when the background orbifold space is static, and 
the ^-derivative term (the third term), which gives graviton's mass as we shall see below. 
So, these two terms reflect the effect of the dynamics of the background branes model on 
gravitational waves. 

Provided that H T (t,y) = T(t)Y(y), the equations of motion (|TTD are reduced to the 
following set of equations for Y(y) and T(t): 



m x 2 



>f +(-) C(y) 3 Y l = 0, (13) 

\CLo/ 

fi + (2A{t) + B{t)) fi + k 2 T x + mV^W^^T, = 0, (14) 
where the dimensionless constant m is defined as follows: 

m =\—^y (15) 

Here, I represents the level of inhomogeneity in the orbifold direction, and m takes discrete 
values labeled by an integer n, as we shall see in the next subsection. Note that the signature 
of the second term in LHS of Eq.(^) is chosen so that the solutions of Eq.(^) satisfy the 
boundary condition (|12j). 



In the case of B = 0, the last term in Eq.(14) provides the eigenvalue (times e /oq) °f 



the dAlembertian on the four- dimensional brane. Therefore the m = modes behave as the 
massless mode in the brane when the orbifold space is static. In this sense, we shall refer to 
m — (I — 0) and m ^ (i ^ 0) modes as "massless" and "massive" modes, respectively. 
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B. Solutions of y-dependent part 

In this subsection, we shall give the solutions of Eq.(J13|) which satisfy the boundary 



condition ([12]). 

For m = modes, Eq.flTB]) reads 



Y l (y)=C 1 + C 2 y, 



(16) 



where C\ and C 2 are integration constants. From the boundary condition (0), we find that 
C 2 must vanish. Therefore, the zero-mode(/ = 0) solution for ^/-dependent part is 

Yi(y) = Cx. (17) 
On the other hand, for m^fl modes, the solutions of Eq. ( |T3"D for y > are given by 



Y l {y)=C{y) 1 ' 2 



C^r-Ciy) 



21 



5/2 



C 2 H[%r-C(y) 



21 



5/2 



(18) 



where C\ and C 2 are constants. 

From the boundary condition at y = 0, we find that the ratio of C\ and C 2 becomes 



9i 
c x 



(19) 



for each I. From the boundary condition at y — np, we obtain 



rr(l) 
-"-4/5 



H (i) 

-"-4/5 



z=2Z/5 



(20) 



=2«( v / 2avrp/3+l) 5 / 2 /5 



This gives the value of I for each excited mode. 
For the case of 21/5 S> 1 , I can be written as 

5 



n it 



(V2a-np/3 + 1 



5/2 



- 1 



(21) 



where n — 1, 2, • • • . airp <C 1 is satisfied in the context of five-dimensional Hofava-Witten 
theory, except for inflationary epoch ||. In this case, Eq . fl2"ID reads 

ao n 



m 



P 



Thus the reduced KK "mass" defined by 

rriKKJ = 



m 



doe 



B(t) 



C(y) 



(22) 



(23) 



takes a typical value of the orbifold energy scale. 
The norm squared of Yi(y) is then given by 



\Yi{y)\ 



2 ^ 2 

w cos 



f \{V2a\y\/3 



5/2 



(24) 



Therefore |Yj(|/)| 2 has peak at each boundary (at y — and y = irp). 



6 



C. Solutions of i-dependent parts 



In this subsection, we shall examine the behavior of solutions to Eg. QT3j). For simplicity, 
we omit the suffix I hereafter. 

For the m = modes, Eq.(|l4]) is exactly solved to yield 



T = D 1 H i u\kt) + D 2 H [ u\kt) 



(25) 



where D\ and D 2 are constants. In contrast, for the m^fl modes, Eg.(|T4"D cannot be solved 
exactly. Therefore, we analyze the behavior of solutions by means of the WKB method by 
rewriting Eg.flI4|) as 



(t 1 / 2 T)-+S(t) 2 t 1 / 2 T = 0, 



(26) 



where S(t) is defined by 



S(t) 



'k 2 



1 

It 2 " 



m 2 t 



2+1-35 



(27) 



In the region where |5| <C S 2 holds, we can use the WKB method to obtain the approximate 
solution 



T{t) ^D 1 (A;)(t5(t))- 1/2 exp 
Since S / S 2 is written as 

\S\ 

s 2 



S(t')dt' 



D 2 (k)(tS(t))-^ 2 



exp 



S{t')d£ 



(28) 



-l + 2(l-35)m 2 t 3 ( 1 - 5 ) 



;i + 4(H) 2 + 4m 2 t 3 ( 1 - 5 )) 3/2 ' 



(29) 



the WKB approximation is good in the region where t m -2 /^ 1-5 )) or t ^> I /(2k). The 
former relation, which is equivalent to vh^yJ S> 1, is satisfied when the time scale is larger 
than the orbifold radius, while the latter is satisfied when the wavelength is shorter than 
the Horizon radius, i.e., k 3> kjj. 

In the (|) background case, since t 1_3<5 is a monotonically decreasing function, the mass 
term becomes negligible compared with the /c 2 -term on the right hand of Eq. (|27j) . Therefore, 
solutions to Eq.flSBD are well approximated by the solutions Eq.(^) in the massless case in 
a sufficient late time for any fixed k. In contrast, in the (j) background case, m 2 t 1-35 term 
increase with time, and the solutions deviate from those for the massless case in late times. 
In particular, for 4m 2 t 3( - 1 ~ <5 ' ^> 1 + 4(H) 2 , the WKB solution is given by 



T^imt^-^y^lD^k) 



exp 



— i- 



2mf 3(l-5)/2 

3(1-5) 



+ D 2 (k) exp 



2mt 3(l-5)/2 

3(1-5) 



(30) 



after an appropriate redefinition of the constants D\ and D 2 . Note that the late time solution 
(1301) does not depend on the wavenumber k explicitly, and the argument of the exponential 
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function in (|30|) is proportional to ttikkT. In particular, from Eq.(|25|) and Eq.(^) or from 
Eq.(p7|) and Eq. (^8|) , we see that the ratio of the amplitude of a massive mode to that of a 
massless mode behaves as 

^ oc S{t)- 1 ' 2 ~ e ^- A ^' 2 = t^> A (31) 
T m= o 

for rriKK ^> k/a. This apparently shows that the massive modes become negligible in late 
times. However, if we consider their backreaction, the conclusion changes significantly, as 
we will see in the next section. 



IV. BACKREACTION OF THE PERTURBATION 

In this section, we study the back reaction of the tensor perturbation on the bulk back- 
ground geometry and on the intrinsic geometry of the brane with the help of the second-order 
perturbation theory. 

First note that if we expand the deviation of the bulk geometry from the background in 
terms of some small parameter, the second-order part 5 2 g satisfies the equation 

(£ {1) 8 2 g) M N = 4(T GW mn + 5 2 T MN ), (32) 

where is the differential operator for the metric perturbation obtained from the linear 
perturbation of the Einstein equations, T GW mn is the effective energy-momentum tensor 
for the linear perturbation 5\g of geometry, which is quadratic in 5ig, and S 2 T MN is the 
second-order perturbation with respect to the field V of the bulk energy-momentum tensor 

4Tmn = ^d M Vd N V - \g MN (^(W) 2 + . (33) 

The explicit expression for T gw mn is given in the Appendix. In contrast to the linear 
perturbation, the spatial average in the 3-dimensional sense does not vanish in general and 
produces a spatially homogeneous contribution to 5 2 g. This contribution can be regarded 
as the backreaction of the perturbation on the background geometry. 

In particular, for the tensor perturbation in the models considered in the present paper, 
the effective energy density pew is given by 

2/gIpgw = -24(T GW °o) (34) 
= {ale 2A ^C{y))-\\H T \ 2 + k 2 \H T \ 2 ) + e~ 2B \H' T \ 2 + • • • , (35) 

under the normalization of the tensor harmonics as (Ty-TP- 7 ) = 1. The leading term of pew 

is given by the first term, which behaves as e~ 2A WHj> ~ t~ 2+s S. In the meanwhile, the 
leading term for the energy density of the V field determining the bulk background geometry 
is given by the potential energy a 2 / (SV 2 ) oc t~ 2S , and its second-order perturbation is given 
by 

n 2 6 T° - 2 62V n 2 T° + & ^ + * ^ ~ 2^ n 2 T° (36) 
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Comparing these energy densities, we find that the effect of pew becomes negligible in the 
(|) background case if S 2 V/V is small. Since the field equation for V is given by 

2a 2 

U gHS2g) (ln(V + (5 2 V))) + 3(V+(W)2 = °> (37) 

up to the second order, {82V) /V is small if fag) is negligible. Hence, in this case the 
backreaction can be neglected in late times. In contrast, for the (f)-model, the decrease of 
Pgw is slower than — Tq . Hence this model is unstable against the backreaction of the tensor 
perturbation. 

Next, we consider the backreaction effect on the intrinsic geometry of the brane. Since 
a full treatment of this problem is very difficult, we only make rough estimate using the 
Hamiltonian constraint along the brane, 

= -2k\T ±1 _ + K 2 - K^K^, (38) 

where is the Ricci scalar of the four-dimensional metric g^ u of the brane, T ±± is the 
component of the energy-momentum tensor along the unit normal to the brane, and K^ v 
is the extrinsic curvature of the brane. As explained above, if we take into account the 
backreaction of the tensor perturbation on the brane geometry, Tj_j_ should be replaced by 
T| + T GW 5 + (^(TljJ in the second-order perturbation framework. Here, K£ is related to 
the intrinsic energy-momentum tensor ^^T^ of the brane by the junction condition 

K = \4 - = T^K- (39) 

Further, the boundary condition on V at the brane is expressed as 

V±V = ±V2a. (40) 

Hence we obtain 

(4 >fl(9 + <**» = (t) R(9) - 24<T° W i} + ^ " ^) , (41) 

where from the Appendix (T GW l) is given by 

4(T™1) = -JL (Hi _ k 2 H^j + GlH 2 + A(i + B)H T H T - ^H T H», (42) 

with b = e B ^C 2 and 

q5 _ J_ (V? 1 ^ _ ^ _ ^_ s J^ s _ 3 

^ 4b 2 V 2 + 4a 2 V 2 QV 2 3C 6 + 4a 2 C ' 1 ' 

By putting the asymptotic estimates for Ht into this expression, we find that in the (J.) 
background case, the H^/a 2 term dominates and decays as t~( 3+s ^( 3 ~ s ^ = t~ 1+2 I <5 I/( 3+ I <5 D. 
Since the background value of ^R(g) decreases in proportion to 1/r 2 ~ t 5 ' 3 , (T GW5 5 ) de- 
creases more slowly than ^R(g). Thus in this case (T GW |) dominates the background value 
for ^R(g) in the late stage, and in order for the FRW nature of the brane to be preserved, 
8 2 V/V must become much larger than unity. This implies that the backreaction of the tensor 
perturbation significantly modifies the evolutionary behavior of the 4-dimensional universe 
on the brane. 
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V. CONCLUSION AND DISCUSSION 



In the present paper we have studied the evolution of gravitational wave perturbations 
in the dynamical FRW brane-world models of the five-dimensional Hofava-Witten theory. 
As the background spacetime, we have used two cosmological solutions, i.e. the (j) and 
(|) solutions, in which the branes represent the spatially flat expanding FRW universes. 
The most important feature of these solutions was the fact that we can solve the evolution 
equation for the tensor perturbation with help of the method of separation of variables in 
spite of the dynamical nature of the brane. Thus we were able to study the spatial and the 
temporal behavior of the tensor perturbation separately. 

Since the model is compact in the fifth dimension, we obtained a discrete spectrum for 
the separation constant which can be interpreted as the graviton mass, and wavefunctions 
for the massless modes and for the massive modes were decoupled as in the case of the static 
brane solutions. However, we found that the spatial behavior of the wavefunctions for the 
massive modes is different from that in the Randall - Sundrum model |13] : in our case they 
have maximums at the branes, but in RS model they have minimums at the branes when 
expressed in terms of the variable Ht adopted in the present paper. This suggests that the 
coupling between the massless modes and the massive modes on the branes may become 
important when we consider non-linear corrections in the models considered in the present 
paper. 

Although the spatial behavior of perturbations for the two solutions was exactly the 
same, their temporal behavior was quite different. Namely, we have found that in the (f ) 
background the temporal behavior of massive modes approach that of massless modes in 
a late time, while in the (j) background the massive modes decay more rapidly than the 
massless modes. We can understand this difference as being caused by the difference in 
the behavior of exp[A(t) — B(t)}, i.e., the difference between the expansion rate of the 
four- dimensional FRW brane and that of the orbifold space. Roughly speaking, waves 
become massless when Xitlkk *C 1, while they become massive when Xtukk ^ 1, where 
A = aoe A ^ \JC(y)k~ 1 is the reduced proper wavelength. Therefore, in the (y) background, 
since exp[A(t) — B(t)] is a decreasing function and \m KK — > 0, every mode becomes 
effectively massless in the late time. On the other hand, in the case of the (J,) background, 
exp[A(t) — B(t)\ is an increasing function, and the modes becomes more and more massive 
with time. Then, as the WKB approximation shows, they suffer from an extra damping in 
proportion to l/(\m KK ). 

This result shows that in both models the tensor perturbation is dominated by massless 
modes. In models such as the Randall- Sundrum model in which the bulk geometry is deter- 
mined by a cosmological constant, massless modes of the tensor perturbation are expected 
to have no important effect on the bulk geometry. In contrast, in the (f) solution of the 
Hoava-Witten theory, the energy density of the bulk spacetime decreases in time. Hence, 
the backreaction of the energy density of the tensor perturbation may become important. 
In fact, we have shown that in the second-order perturbation framework, the contribution of 
the tensor perturbation supersedes the original background energy density determining the 
bulk geometry in the (f ) case. Hence this background solution is unstable against non-linear 
corrections. Although this result was obtained for a special brane motion obtained under 
the assumption that the brane contains no matter apart from the <fi field, it also holds for a 
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more realistic brane which contains ordinary matter. It is because the essential feature of the 
temporal behavior of the tensor perturbation does not depend on the boundary condition 
at the brane. 

We have also examined the 2nd-order backreaction of the tensor perturbation on the 
intrinsic geometry of the branes, and have found that the backreaction effect significantly 
alters the time evolution of the brane geometry for the (j) solution, although the backreac- 
tion on the bulk geometry is negligible for this solution. This result is consistent with the 
naive expectation that the massive modes of the tensor perturbation behave as dark matter. 

These results suggest that the stability against the backreaction can be used as a crite- 
rion to physically acceptable brane-world models and to discuss cosmological implications 
of models. Thus it will be interesting to analyze the non-linear stability of the Randall- 
Sundrum models as well as of more realistic solutions in the Horava-Witten theory. 
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APPENDIX: 



In this appendix, we give the expression for the second order part with the respect to 
the perturbations of the (m + n)-dimensional Einstein-Hilbert action 



S=— I d m+n x^(R-2A), 



(Al) 



and its energy-momentum tensor, where k and R are the (m + n)-dimensional gravitational 
constant and Ricci scalar, respectively, and A is the cosmological constant. 

By decomposing the metric g pv into the background g^ u and the perturbation h pu as 



9tiv = g^u + h 



(A2) 



and substituting it into Eq.([Al|), we obtain the following expression for the second order 
part with the respect to the perturbation h pv of the (m + n)-dimensional Einstein-Hilbert 
action: 



im+n 



+ -(h 2 - 2h^h^)R + -(2W p h p u - hh^)R pv - -A(h 2 - 2h^h 



(A3) 



By taking the variation of this action with respect to the background metric g pu , 

5S 2 = J d m+n x^^5g^, (A4) 
we obtain the following energy- momentum tensor T pu \ 
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Ak 2 T, 



pa 



1 



2 v- — - W 



(/i 2 - 



+ h P °; U h L 



Oh P\~\h - 2h pa h - 2h pa h + 2/? p /? CT 4- 2h p h a 



\p lb pa:v 1 — "''p/i \pp ib vu \^ ,v V(T\p 1 ^ J,t ' fi;va" j p \ p\v ib p\a 

+4h pa h 

+QRp p axK x h% + 2{h p i K v - h^h p<7 ) ;pa -(R- 2A)(/iV " 2^p^ p ,) 
_ 2h pa h pa ), yi/ - AR^h/h/ - hh p v ) - 2R pa (h p v h p ° - h pa h p 



"jXV ) 



(A5) 



where □ is the (m + n)-dimensional d'Alembertian. 

Under the notation adopted in JHJ, the unperturbed background geometry in brane- 
world models is expressed as 



ds 2 = g^dz tM dz u = g a b(y)dy a dy b + r 2 (y)da, 



(A6) 



where the metric 



da 2 = r )ij{x)dx l dxi 



(A7) 



is that of the n-dimensional space with a constant sectional curvature K. The tensor mode 
of the metric perturbation is expanded as 



h 



ah 



0? hai 0? hij 



2r 2 H Tij . 



(A8) 



For this tensor perturbation, the spatial average of the energy-momentum tensor 
the spatially flat (K = 0) case is given by 



in 



ab 



k 2 

{DH T ) 2 + —Hi \ + -H T Dr ■ DH T + AH% 



4 



-D a H T D b H T - 2H T D a D h H T - -H T D a rD b H T + G ab H 2 



K 2 {T GW )=r% 



^{(DH T ) 2 + ^H 2 } + (A + -Gt)H, 



- 2{r 2 (DH T ) 2 + k 2 Hl}+Ar 2 ^G k k + AjH, 
where we have normalized as (TTyT^) = 1 and (T^T j) = Cij. 



(A9) 
(A10) 

(All) 
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